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Abstract 

Context. A number of microlensing dark-matter surveys have produced tens of millions of light curves of individual background 
stars. These data provide an unprecedented opportunity for systematic studies of whole classes of variable stars and their host 
galaxies. 

Aims. We aim to use the EROS-2 survey of the Magellanic Clouds to detect and study the population of beat Cepheids (BCs) in 
both Clouds. BCs pulsating simultaneously in the first overtone and fundamental modes (FO /F) or in the second and first overtone 
modes (SO/FO) are of particular interest. 

Methods. Using special software designed to search for periodic variables, we have scanned the EROS-2 data base for variables 
in the typical period range of Cepheids. Metallicities of FO/F objects were then calculated from linear nonadiabatic convective 
stellar models. 

Results. We identify 74 FO/F BCs in the LMC and 41 in the SMC, and 173 and 129 SO/FO pulsators in the LMC and SMC, 
respectively; 185 of these stars are new discoveries. For nearly all the FO /F objects we determine minimum, mean, and maximum 
values of the metallicity. 

Conclusions. The EROS data have expanded the samples of known BCs in the LMC by 31%, in the SMC by 110%. The FO/F 
objects provide independent measures of metallicities in these galaxies. The mean value of metallicity is 0.0045 in the LMC and 
0.0018 in the SMC. 

Key words. Magellanic Clouds - Cepheids - Stars: fundamental parameters 



1. Introduction 

Over the last decade, large microlensing surveys have been 
conducted to search for planetary- to stellar-mass baryonic 
dark matter in the Milky Way. While, interestingly, th e 
original search was in fact negative (Tissera nd et alj |2007). 
the resulting tens of millions of homogeneous, well-sampled 
light curves of the background stars have stimulated many 
new programmes on variable stars. In particular, specific 
stellar populations like Cepheids can now be studied far 
more thoroughly, particularly in the Magellanic Clouds. 
Among these objects, the subsample of double-mode or 
beat Cepheids (hereafter BCs) is of particular importance 
in gaining information on the structure and evolution of 
massive stars. 

According to th e General Catalogue of Variable Stars 
(jSamus et al.l [2003 ). the Galactic sample of known BCs 
is rather small: up to now, only 16 objects of this type 
have been identified. In the Large (LMC) and Small (SMC) 
Magellanic Clouds the situation is more favourable, thanks 



in pa r ticula r to the wor k of the EROS-1 (fBeaulieu et al.l 
119951 119971b MACHO (lAlcock et all 119951. Il999h . and 



* Based on observations made by the EROS-2 collabora- 
tion with the MARLY telescope at the European Southern 
Observatory, La Silla, Chile, 
deceased 



OGLE lUdalski et all 119991 : ISoszvnski et all I2000T) groups. 
EROS-1 detected one BC in the LMC and 11 in the SMC, 
while the MACHO group found 45 BCs in the LMC. The 
OGLE collaboration reported the discovery of 93 BCs in 
the SMC and 76 in the LMC, very recently extended by 
the OGLE-III data t o a region of 40 deg 2 of the LMC 
(jSoszvnski et al.ll2008l ). This thorough analysis yielded 266 
BCs of different types discussed below. The present study 
will report results on an even larger region of 88 deg 2 in 
the LMC and 10 deg 2 in the SMC. 

A BC pulsates either in the first overtone and funda- 
mental modes (FO/F), or in the second and first overtone 
modes (SO/FO jB The studies listed above clearly estab- 
lished that the period ratio (higher to lower mode) of the 
FO/F pulsators is around 0.72, while that of SO/FO objects 
is near 0.80. For the FO/F pulsators it is well known that 
the period ratio depends on the metalli city Z, as seen e.g. 
in the Petersen diagram (|Petersenlll973f) of the period ratio 
vs. the period of the lower mode. From OGLE photometry 



Two c ases are known of pu lsations in the first and third 
overtones (|Soszvnski et al.ll2008l ). 
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and stellar atmosphere models. iKovacsl (|2008T ) showed that, 
in both Clouds, the average metallicity of FO Cepheids is 
lower than those pulsating in the fundamental mode. 

The chemical enrichment history of t he Magellanic 
Clouds has been recently reexamined by ICarrera et al.l 
( 2008a. b), who used the infrared Call tripl et as a metal- 
licity indicator for red g i ant b ranch stars (|Carrera et al.l 
l2007l) . [Leisv fc Dennef eld (2006) used planetary nebulae to 
trace the heavy elements in the Magellanic Clouds. They 
concluded that only Ar is a suitable tracer of the metal- 
licity of the progenitor stars, instead of a n average of sev- 
eral different elements. Ildiart et al.l (|2007f ) also studied the 
chemical evolution of the SMC from planetary nebulae and 
concluded that the star formation histories of the LMC and 
SMC are different from each other and in many respects still 
controversial. 

The situation is complicated by the fact that metallicity 
measurements in the Magellanic Clouds yield many differ- 
ent values, depending both on the location and age of the 
stars studied. Thus, a single firm value of the metallicity in 
these galaxies cannot be assigned. Abundances of a-, iron- 
peak and s-process elements in individual stars are now 
being me asured with high-resol ution spectrometers, for ex- 
ample by iPompeia et alj ([2008), but it remains important 
to contribute to these efforts by adding a new and indepen- 
dent way to measure th e metallicity, as we do below. 

In a similar study, iBeaulieu et alj (|2006f ) deduced a 
metallicity gradient in M33, only from the pulsational prop- 
erties of 5 newly discovered BCs in M33 and with the help 
of stellar pulsation theory and mass-luminosity relations 
derived from evolutionary tracks. The result was in good 
agreement with the standard spectroscopic metallicity gra- 
dients as determined from H II regions, early-B supergiant 
stars, and planetary nebulae. BCs are thus a new and in- 
dependent probe of galactic metallicity. In this paper we 
present the results of a systematic search for and analysis of 
BCs in the Magellanic Clouds from the EROS-2 database. 
Numerous previously unknown objects are detected, and 
metallicities are derived from the properties of FO /F BCs 
in both Clouds. 

2. Observational setup and BC search 

The EROS-2 experiment was conducted between July 1996 
and March 2003 using the dedicated MARLY telescope 
(lm Ritchey-Chretien, f/5.14) at ESO, La Silla. A dichroic 
beam-splitter allowed simultaneous imaging in two broad 
non-standard passbands, £?e (4200-7200 A, "blue" chan- 
nel) and R E (6200-9200 A, "red" channel). Each camera 
contained a mosaic of eight 2K x 2K LORAL CCDs with a 
pixel size of 0.6"and a field of view of 0.7°(right ascension) 
x 1.4° (declination). 

The light curves of individual stars were constructed 
from fixed positions on templates using special-purpose 
software for p hotometry of EROS-2 images called PEIDA 
(|Ansaril IT996) . This resulted in a database of tens of mil- 
lions of light curves, which will be publicly available in the 
future from the CDS through a dedicated server. For a typ- 
ical object, 800 to 1000 epochs of observation are available. 

Each light curve is named according to the rules given 
bv iDerue etaTI (|2002t ). As an example, the star J053929- 
703819 with J2000 equatorial coordinates of 05:39:29.15, 
-70:38:19.2 is also referred to by the internal EROS-2 num- 
ber lm059115619 (the 5619 th star in the field lm059, CCD 



1, quarter 1). In the following, we refer to the objects by 
their internal number. 

A systematic search for periodic objects using the -Be 
observations was performed among stars detected as poten- 
tially variable by an automatic time series analysis p i peline , 
base d on software developed by IBeaulieu et alj (|l997l ) 
and ISchwarzenberg-Czernvl (|2003l ) , and on an Analysis of 
Variance (AoV) method. For each periodic pulsator we con- 
structed a Fourier decomposition of the light curve, i.e. the 
signal wa s represented by R( t) = Rq + Y^iZi Ri cos(^i(t — 
t ) + <Z> z ) (jSimon fc Ledll98lh . The quantities R kl = R k /Ri 
and <3?/ci = — k <I>i (1 < k < 3) were defined as usual. 

The maximum of t he st atistic ®AoV as defined by 
ISchwarzenberg-Czernvl ()2003l ) was used as a quality pa- 
rameter for the periodograms, which were stored when 
&AoV ^ 20. In this way, a total of about 80 000 variables 
was found in both Clouds. The Fourier model was then sub- 
tracted from the observational data before a second step of 
period search by the same pipeline. Double-mode objects 
thus emerged quite clearly from the resulting periodograms. 
Potential BCs were identified by having period ratios near 
0.73 and 0.80 for FO/F and SO/FO pulsators, respectively. 

To construct the final samples, the light curves of all 
candidates (folded according to the first mode of pulsation) 
were inspected visually, using the versatile tool TOPCAT 
developed by the Starlink consortium^ and customized to 
display SuperMongo plots. On a typical BC light curve, 
folded according to either the first or second period, the 
beat phenomenon is easily recognised because the light 
curve is broader than that of a single mode Cepheid. This 
was our ultimate criterion of selection for objects with the 
lowest values of OaoV- 

Figure [1] shows six typical LMC objects, alternately of 
FO/F and SO/FO type. The two first are firm detections 
with high values of QaoV for the second pulsation mode, 
the two in the middle are detected more marginally, and 
the two last have been rejected as having too low values 
of QaoV- The left-hand panels show the folded light curves 
(blue magnitude vs. period of the second mode of pulsation, 
in days), after subtraction of the Fourier model of the first 
mode. The right-hand panels show the corresponding peri- 
odograms (&AoV vs. frequency of the second mode period). 
In each of these periodograms the highest peak corresponds 
to the most probable frequency. 

3. General analysis and determination of the 
metallicities 

The visual inspection yielded 74 FO/F pulsators in the 
LMC and 41 in the SMC as well as 173 SO/FO stars 
in the LMC and 129 in the SMC. The pertinent data 
for these stars are given in separate tables for the FO/F 
and SO/FO pulsators in each galaxy (LMC and SMC) 
which are only available in electronic form at the CDS via 
anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-l-A/. Table Q] 
summarizes their contents. 

Figure displays the spatial distribution of our BCs 
in both Clouds. The bar of the LMC is clearly visible, al- 
though FO/F and SO/FO objects are widely distributed 
in both the LMC and SMC. We note a surdensity of 
Cepheids in the northern part of the LMC (75° < a < 90°, 

2 Available at http://www.star.bris.ac.uk/~mbt/topcat/ 
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Figure 1: Folded light curves and corresponding pcriodograms of 6 typical LMC BCs, alternately FO/F and SO/FO 
objects, from the most significant (top) to non-significant secondary modes (bottom). Left: Be magnitude vs. period 
(days); right: QaoV for the second pulsation vs. frequency (days -1 ). The name of each object is given in the right-hand 
panels, together with the maximum QaoV for the first pulsation mode. See text for details. 
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Table 1: Contents of the data tables for the new BCs in the 
LMC and SMC. 
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6 for FO/F pulsators only 

c superscript i refers to the i th pulsation mode 



—68° < S < —66°), corresponding to a region of star for- 
mation. We have checked that BCs have a spatial distribu- 
tion similar to that of single-mode Cepheids by performing 
Kolmogorov-Smirnov (KS) tests using the stats package 
of R, a language and environment for statistical computing 
(|R Development Core T eam 200j|) . We computed the mean 
coordinates of the barycentre for the subsamples of F and 
FO objects, obtained to first approximation from a visual 
cut in the period-luminosity relation, and for the BCs. We 
then calculated the distance to the barycentre for each ob- 
ject and made the one-sided KS tests reported in Tabled 
on the resulting distance distributions. 

We first verified that the distributions of F and FO 
Cepheids are almost identical (first and third lines in Table 
[2]). Note however the relatively high p- value (i.e. the prob- 
ability to reject the null hypothesis) of this KS test in the 
LMC. The low p- values in the second and fourth lines in 
Tabled show that BCs and single mode Cepheids are dis- 
tributed similarly. The D statistic represents the maximum 
distance between the distance distributions of populations 
(a) and (b). 

We have compared our findings with the most recent 
studies from the OGLE group. Table [3] summarizes the 
numbers of stars which are in common between EROS-2 
and OGLE, and those which were found only by EROS- 
2 or OGLE, respectively. For both the FO/F or SO/FO 
samples, the first column list the number of objects de- 
tected by both surveys, while the next two columns give 
the number of stars t hat are unique to each . OGLE data 
for the LMC are fro mlSoszvnski et all (|2008h . for the SMC 
from lUdalski et all (|1999T) . Although a detailed analysis of 
this comparison is beyond the scope of the present pa- 
per, it appears that a large fraction of BCs are detected 
by both experiments, especially in the L MC for which an 
extended OGLE-III sample is available (jSoszvnski et al.l 



2008). The present work has, however, significantly in- 
creased the known samples of BCs, by some 31% in the 
LMC and 110% in the SMC. 

There can be many reasons why some objects appear in 
the EROS-2 data and not in OGLE and vice versa. E.g., 
the spatial coverage is not the same, and the EROS-2 pho- 
tometry is PSF-based while OGLE uses a differential image 
analysis method. In addition, the periodicity search algo- 
rithms are different and it should be useful to study whether 
they are sensitive to different numerical effects or not. In 
the LMC, 66 (22 FO/F and 42 SO/FO) of the 82 objects 
that are detected by EROS-2 only are located in regions 
out of the OGLE-III fields. In the SMC, 69 (17 FO/F and 
52 SO/FO) of the 81 objects identified only by EROS-2 are 
external to the OGLE-II fields. Inversely, on the 127 BCs 
found by OGLE only, we estimated than 40% are located 
either outside or on the edge of the EROS-2 CCDs, which 
did not allow a proper photometric follow-up. For the 60% 
remaining, a mixed between EROS-2 lower photometric ac- 
curacy and difference in algorithm sensitivity is thought to 
be the answer. Clearly, this question should be examined 
in more details in the future. 

Figure [3] displays the period ratio Piq = P\/Pq vs. 
logPo, a Petersen diagram with all FO/F and SO/FO 
EROS-2 objects. While SO/FO Magellanic objects are 
mixed in the upper part of the diagram, a clear separation 
between FO/F LMC and SMC BCs is visible in th e lower 
part. This was previously noted by iBeaulieu et al.l (|l997f ) 
who explained from a theoretical point of view that the 
observed period ratio results from a combination of two 
opposite effects: (i) the lower the metallicity, the smaller 
the opacity bump which affects the pulsation; (ii) at lower 
metallicity the structure of the star changes: at fixed mass, 
the luminosity increases and, therefore, the period ratio de- 
creases. Moreover, the position of the blue loop of an evo- 
lutionary track of a Cepheid with respect to the instability 
strip is different in the LMC and in the SMC because of 
their different metallicity. As a consequence, lower pulsa- 
tion periods are favoured in the SMC, while the highest 
periods are found in the LMC. 

The m etallicity of FO/F BCs has been recently reex- 
amined bv lBuchler fc Szabol (|2007f ). It has been known for 
some time that the location of a BC in the Petersen di- 
agram correlates strongly with metallicity Z. We plotted 
all our LMC and SMC BCs on such a Petersen diagram in 
Figs, g] and H 

iBuchler fc Szabol (|2007t ) developed a convenient tool for 
determining metallicity limits, and Z max from the 

location of a BC in the Petersen diagram, based only on 
accurate observed pulsation periods. From linear nonadia- 
batic convective stellar models they constructed a compos- 
ite Petersen diagram (Figs. [4] and [5]), in which the isometal- 
licity lines delimit the lower and upper boundaries, respec- 
tively, of the range in which both F and FO modes are 
linearly unstable; this is the linear criterion for beat pulsa- 
tions. 

Thus, if a given BC falls below the upper period bound- 
ary for a certain value of Z in Fig. 21 and above the cor- 
responding lower period boundary for the same Z in Fig. 
then this is a possible value of Z for this star. Thus, 
for example, for the LMC BC with log P = 0.521 and 
Pxo = 0.7067, Fig. H yields Z mm - 0.0077, while Fig. 
gives Z max ~ 0.0093. It is interesting to note that the ob- 
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Figure 2: Spatial distribution of BCs in the Magellanic Clouds. Open circles: FO/F objects; filled circles: SO/FO objects. 



Table 2: Results of the one-sided Kolmogorov-Smirnov tests performed on the spatial distributions of two different sets 
of Cepheids, noted (a) and (b). See text for further detail. 





Set (a) 


Set (b) 


D b 


p-value c 


LMC 


F (1772)" 
F + FO (2679) 


FO (907) 
FO/F + SO/FO (247) 


0.0455 
0.0994 


0.1660 
0.0245 


SMC 


F (1256) 
F + FO (2537) 


FO (1281) 
FO/F + SO/FO (161) 


0.0568 
0.1322 


0.0335 
0.0107 



"Numbers in parentheses are the sizes of the populations 

b D defines the maximum distance between the distance distributions of (a) and (b) 
c The p- values are the probability that the distributions are different 
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Table 3: Comparison of the EROS-2 (this work) and OGLE BC samples in the LMC and SMC. 





FO/F 


SO/FO 




EROS-2 & OCLE 


EROS-2 


OGLE 
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Figure 3: Petersen diagram of FO/F EROS-2 BCs (lower; 
filled dots: LMC; crosses: SMC) and SO/FO objects (upper; 
open circles: LMC; stars: SMC). 



servational data follow the theoretical isometallicity curves 
rather than a simple linear correlation. 

In practice, one interpolates among the isometallicity 
lines with a simple code. Using this approach, we have 
computed the metallicities of our BCs, and give the mean 
Z-mean together with Z m i n and Z max in the electronic ta- 
bles. Figure [6] shows the distribution of mean metallicity for 
the FO/F BCs in the Magellanic Clouds. The size of the 
bins is 3.25xl0 -4 , corresponding to the mean uncertainties 



in Z„ 



which is 5.2 x 1CT 4 in the LMC and 2.3 x KT 4 in 



the SMC. 

A comparison of these data with other determina- 
tions is not straightforward, becau se previous spect ro- 
scopic measurements su c h as t hose by ICole et all (|2005l ) or 
ICavrel de Strobel et all (|200Lh report [Fe/H] as obtained 
for stars from different populations, thus from different 
stages in the chemical evolution of the host galaxy than 
that of the young BCs. Also, we determine Z 1 not [Fe/H], 
and allowance f or the non-Solar ab undance ratios in the 
LMC and SMC (jPompeia et al.ll2008l ) would be required to 
make a direct comparison. 

Three low-luminosity objects with periods ~ 0.5 d fall at 
the left boundary of the Petersen diagrams: Im0316m27403, 
Im0470n7559 and Im0682ml5849 (Fig. [3]). They also appear 
to have Z > 0.026, which is not realistic. This therefore 
seems to rule out that they are very low mass Cepheids i n 
the first crossing phase, as argued bv lBaraffe et al.l (|l998h . 
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Figure 4: Composite Petersen diagram (Pro = Pi/Po vs. 
P Q ) of the FO/F BCs from EROS-2 compared to other 
known stars. The lines delimit the periods for which both 
F and FO are linearly unstable at given Z (lower bound- 
ary in Z for given P 10 and P ). The metallicity increases 
downward from Z — 0.001 (top line) to 0.014 in steps of 
0.001. Red rhombuses: EROS-2, LMC; blue squares: EROS - 
2, SMC; black stars: Gala ctic objects (ISamus et alj|2003 ): 
pink circles: M33 objects (|Beaulieu et al J 120061 ). 



In fact, they may not be bona fide Cepheids within the 
range of our calculated models. However, we recall that the 
only assumption in our computations is that they are cen- 
trally condensed, like Cepheids, so that they can be repre- 
sented by envelope models. It could therefore be that these 
stars are less evolved. Concerning the 4 BCs (2 LMC + 2 
SMC) not represented on Figs. 0] and [5] as they are outside 
the range of Pro on that figures, we did not find any partic- 
ular behaviour. Therefore we suggest, that these BCs may 
be used to localized zone in the Magellanic Clouds of lower 
and higher metallicity than average. 
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Figure 5: As Fig. [4] (same symbols), but showing the upper 
boundary in Z for which models are unstable for given Pq 
and Piq. 



n 



Figure 6: Distribution of mean metallicities of FO/F BCs in 
the Magellanic Clouds. Solid line: LMC; dashed line: SMC. 
The bin size is 3.25 xlO -4 , corresponding to the mean error 

*^f Z mean . 



the nonlinear range is closer to the upper linear boundary 
(Fig. [5]) than the lower one (Fig. |4|). On average, the non- 
linear values of Z are therefore expected to be somewhat 
higher than the linear ones. 



4. Discussion of the uncertainties in Z 



4.2. The chemical composition 



We comment here on the three main sources of uncertainty 
in the inferred metallicities. 



4.1. The use of a linear Beat Cepheid criterion 

In Fig. [7] we show the derived Z values as a function of log 
Po . The error bars reflect the consequences of using the lin- 
ear criterion for beat behavior, i.e. simultaneous instability 
in the fundamental and first overtone modes. A comparison 
of Figs. 0] and O which show the upper and lower boundaries 
of Pio as a function of Z , indicates that the uncertainties in 
Z become very large for the longer-period BCs. It is indeed 
well known that the linear BC criterion is far too relaxed, 
so that real BCs occu py only a subregion of the correspond- 
ing p arameter space (|Stellingwer i 1 1974 iBuchler fc Kovacsl 
119861) . 

The obvious remedy is to perform nonlinear BC model 
calculations. This is a very computer intensive and time 
consuming endeavour, but is now under completion (Szabo 
& Buchler 2009, in preparation). The work to date shows 
that the width of the double-mode instability strip is nar- 
rower by about a factor of 5 in log T in a theoretical HR dia- 
gram. In a diagram like Fig.[5]this translates into a spread in 
Pio of less than 0.0005; moreover, this is relatively uniform 
in log Pq, i.e. does not open up towards the long periods 
as with the linear criterion. It is also found that, overall, 



Our models use OPA L opacities computed with the revised 
solar composition of lAsplund et al.l (|2005f ). This causes a 
substantial change in the inferre d metallicity as compared 
to that obtained with th e older Grevesse fc Noeisi {l993) 
solar composition used bv lBuchieiT ( 2008[ h The inferred av- 
erage ('linear') Z in the LMC is reduced from 0.0070 to 
0.0044 and that of the SMC from 0.0028 to 0.0018. This 
is of course a gross overestimate of the true uncertainty 
in metallicity, because any future revision of the compo- 
sition is expected to be s malle r than the difference be - 
tween lAsplund et all (p005l) and[G revesse fc Noeld (|l99i . 
We also note that Magellanic Cl oud stars in general do 
not have Solar abundance ratios (jPompeia et al.ll2008f ) as 
we have assumed in our models, and we do not know how 
much this might change the results. 

4.3. Modelling uncertainties 

IBuch ler (2008|) has shown that the approximation of enve- 
lope models breaks down for the short-period Cepheids, and 
that it is necessary to compute full stellar models (which 
we are not equipped to do at this stage). One may there- 
fore wonder whether the lower metallicity at shorter periods 
seen in Fig. [7] is a result of this uncertainty in the models, 
or whether there is a physical or evolutionar y reason for 
the lower metallicity, e.g. ICordier "eTall (|200l . 
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Figure 7: Plot of the metallicity Z vs. the period of the fundamental mode. Error bars are derived from the lower and 
upper boundaries in Figs. 0] and [5] 



5. Conclusion 

We have reported a new survey of the BC populations in 
the Magellanic Clouds which co mplements the recent work 
of the OGLE group in the SMC (jSoszvnski et al.ll2008[ ) and 
extends it over a larger area in the LMC. Our new, large 
FO/F sample yields a new metallicity distribution for an- 
other stellar population in the LMC and the SMC. Together 
with previous studies of HII regions, planetary nebulae or 
red giants, this contributes to a more complete view of the 
metallicity of the Magellanic Clouds and its links to stellar 
evolution. We deduce from the present data a mean metal- 
licty of 0.0045 in the LMC and 0.0018 in the SMC. 



The demonstration that such BC stars are a useful 
and independent way to determine metallicities has con- 
siderable promise for future studies in the Local Group. 
In M33, 5 BCs were found among a sample of more than 
2000 Cepheids (|Beaulieu et al.ll2006fl . and we have shown 
here that tens of BCs can be detected in the Magellanic 
Clouds, amon g a total population of 2 000-3000 single- mode 
Cepheids (see ISoszvnski et all (|2008l ) and Table [2] above). 
The frequency of BCs in other galaxies of the Local Group 
is likely to be between these values. M31 is the obvious first 
candidate for such a search. 
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